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(54) Method and system for decoding distorted image and symbology data 



(57) Disclosed is a system and method for reading 
and decoding distorted data, e.g. , bar code label or other 
symbology data, by fitting a scanned signal model /(/, 
9), using a Deterministic Expectation -Maximization 
(DEM) algorithm, to the measured or detected data sig- 
nal y(t). The DEM algorithm, which is a novel variation 
of the conventional Expectation-Maximization (EM) al- 
gorithm, enables the system and method according to 
the invention to determine the values of a plurality of 
parameters 9 defining the scanned signal model /(/, 9) 
as it is being fitted to the measured data signal y(f). 



Based on these values, characteristic information such 
as data elements or symbols encoded within the image 
data are generated using conventional signal process- 
ing techniques. The use of the DEM algorithm advanta- 
geously deblurs the signal by reducing the effects of in- 
tersymbol interference. In an alternative embodiment, 
the DEM algorithm is supplemented by a conventional 
zero crossings decoding technique. In this manner, the 
zero crossings technique is used to obtain initial param- 
eter estimates for the DEM algorithm in determining the 
parameters 9 of the scanned signal-model /(/, 9). 
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Description 

Backeround of the Invention 
5 1. Field of the invention 

The invention relates to methods and systems for processing distorted image and symbology data More particularly 

ErZS " t0 meth0dS SyStemS ,0f readiP9 ' deCOdin9 3nd reCOVedn 9 bar Symbols and otJ^ym: 
bology data from representative signals. &ym 
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2. Description of the Related Art 

The use of symbology data is becoming increasingly prevalent in our society For example, a particular bar codina 
arrangement known as the Universal Product Code (UPC) has been established as the industry standard 
food processing and other retail industries. grocery. 

K, a ,? ene ?H '.K 8 ' ? de J abe,S USe ver1ica,, y a,, 9 ned - alternating stripes of two colors, typically light absorbing usually 
black, pnnted -bars" and light scattering (diffusely reflecting) "white" spaces. The label data is encoded in the reS 
w^ths of the aUernat.ng bars and spaces. Such relative width encoding allows label reading to be relatively nsenl ^ 
to geometnc effects such as label size : reading distance and to a certain extent, label orientation ins " e 
n Hof S "Ti Syslems JyP ica,, y involve ™ optical scanner, a pholodeteclor and associated signal processing (e 
L d C ° n ? n 3 ° r S ° ftWare ' SCannef ' WhiCh ° ften inC,udes a Iaser <* othe ' su "able source is 

y ^ 9 h 3 ! ° f ,GnS SyStem) t0 Pr ° ieCt 3 bri9hl SP ° 1 ° f ,j 9 ht ' hav,n 9 a certa - *™ afa prescribed 
d.stance, that ,s moved across the bar code label. The photodetector, positioned approximately coaxial with the light 
source, measures the t.me-varying intensity of light backscattered from the label and converts it to a correspond 
time-varying electrical signal current. »n«««"g 

» « ^ S ^ em l Si , 9na i pfOCeSsin9 ,hen converts ,he time-varying electrical signal, which is an analog waveform to 
a pu.se width modulated digital signal. The digital signal is then transformed to a binary representation of the daTa ? 

thefe3 t^n. nUmberS f deSCrib j n9 ,he Widths ° f ,he al,ernati "9 ^ °' bar code. The .abe. data is interpreted 
therefrom using any one of a number of conventional decoding techniques / 

of ,hI h ^rT," ab , i K y t0 , C T 6 ' abelS depSndS ° n the dimensions ° f 'he p/ojected light spot relative to the size 
till ? >' ' ' S,,y ° f the Pr ° ieC,ed Spo1 re,a,ive 10 ,he bar aflects.'the systems ability to read because 
substantial noise sources are present. Laser power is limited by eye safety considerations and thus provides an upper 
..mrtonavateblesigna^^ 

inr.nH« 'T """I'r" ?!" pr °> eCted spot introduces corresponding blurring in the bar code image. Such blurring 
includes ■ntersymbol interference (also called convolution interference), which is the effect of light backscattered from 
bars adjacent to the bar of interest. 

r~,H CO T nti0 ?' L ead ^ 9 meth ° dS - WhiCh ,ypiCa " y are based on linear ,ilterin 9 and/or feature detection work well in 
eading bar code data rt the blurring is small relative to the size of the bars. In such cases, reading performance is 

ZZTZ^S ,h : h S ' 9na| - l °: n0iSe ra,i °- Yet ' 95 re,alive b,urrin 9 ™ ease s ^ 'abe- moves furtLr Z focus 
Tn JTrl r ' ndependenl of ,he signaMo-noise ratio. Thus, in these cases, increasing the detected 
signal does not necessarily reduce such reading errors. 

Conventional reading methods are limited effectively to cases in which blurring caused by intersymbol interference 
methods d'c "J 0 * me S ?r °' ? ' abe ' ^ 03565 ^ 15 ™ Ch more — ent ona, readTg 

Bc^^Z syTeT ' nCreaSed ,hat in,erSymb °' ,n,erference blu ™9 " a * on the overall reading 

daJ^™ h° H Vai ' ab ' e 3 Sy5,em 3nd mS,h0d ,0r readin 9 dala such as bar code *» °' o'ber symbology 
^cussed above ' ma9e ° PliCa ' ChafaC,er feC ° 9nUion that improve upon the deterministic distortions 

ar. inH^'nV T Xem , a0d mCth0d are d6Sired ,hat fedUCe intefs y"bol interference and other types of blurring that 
are independent of signal-to-noise ratio. a 

Summary of the Invention 

mp.h^r^ a ", d aCCOrdin 9 to the invenJon are as set out in the independent claims, preferred forms of the 

s~; n e d ^cT^ssr c,aims - in pre,erred ,orms ° f the syste - - — * -J « 

The invention is embodied in a system and melhod for reading and decoding distorted data. e.g.. barcode label 
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data, facsimile and other scanned image data having distortions therein. In particular, it is embodied in a system and 
method that decodes data such as bar code label data by fining a scanned signal mode! / (/, 0), using a Deterministic 
Expectation-Maximization (DEM) algorithm, to a measured or detected data signal y(t). The DEM algorithm, which is 
a novel variation of the conventional Expectation-Maximization (EM) algorithm, enables the inventive system and meth- 

5 od to determine the values of a plurality of parameters 0 defining the scanned signal model /(/, 0) as it is being fitted 
to the measured data signal y(t). Based on these values, characteristic information such as data elements or symbols 
encoded within the image data are generated using conventional signal processing techniques. The use of the DEM 
algorithm advantageously deblurs the signal by reducing the effects of intersymbol interference. In an alternative em- 
bodiment, the DEM algorithm is supplemented by a conventional zero crossings decoding technique. In this manner, 

io the zero crossings technique is used to obtain initial parameter estimates for the DEM algorithm in determining the 
parameters 0 of the scanned signal model j(t t 0). 

Brief Description of the Drawings 

15 in the drawings: 

FIG. 1 is a simplified block diagram of a data label reading method in which aspects of the invention find application; 
FIG. 2a is a schematic diagram of a bar code reading system in which aspects of the invention find application; 
FlG.2b is a schematic diagram of a pulse train output from the bar code reading system of FIG.2a; 
20 FIG. 3 is a detailed view of bar code scanner light as it relates to a portion of a bar code label, 

FIG. 4 is a graphical diagram of the measured data signal y(t) and the scanned signal model j(t t 0) fitted thereto 
using a conventional decoding technique; 

FIG. 5 is a graphical diagram of the measured data signal y(t) and the scanned signal model y(r, 0) fitted thereto 
using a conventional zero crossing processing technique; 
2S FIG. 6 is a graphical diagram of the measured data signal y(t) and the scanned signal model y(f, 0) fitted thereto 

according to an embodiment of the invention; and 

FIG. 7 is a graphical diagram of the measured data signal y(t) and the scanned signal model y(t, 0) fitted thereto 
according to an alternative embodiment of the invention. 

30 Detailed Description / 

In the following description similar components are referred to by the same reference numeral for clarity and con- 
sistency among the drawings. 

The invention is described primarily in terms of a particular class of embodiment, namely, bar code readers and 
35 methods of reading barcode labels. However, the invention is not thus limited, as will be appreciated by those skilled 
in the art. 

FIG. 1 illustrates, in simplified block diagram form, the inventive data label reading method 2 described herein. 
The first step 4 of the method is to scan a label using an appropriate scanner. The label is, e.g., a bar code label, 
scanned image data or other symbology data potentially having distortion therein As discussed previously, scanning 
40 the label involves moving a projected spot of light across the label. The light backscattered from the label is indicative 
of the label's pattern of bars and spaces. 

The next step 5 in the method is to measure the light y(t) backscattered from the label using a photodetector that 
generates an electrical signal based on the intensity of the measured light. The backscattered light varies over time 
as the projected light from the scanner moves across the label. 

Once the measured data signal y(t) has been captured, characteristic information in the measured data signal y 
(t) is determined using a decoding step 6. For example, in the case of bar code labels, such determination includes 
decoding the symbology data encoded in the label data as interpreted from the measured data signal y{\). As mentioned 
previously, conventional decoding steps are based on linear filtering or data features within the measured data signal 
Mt). 

However, the inventive method described herein decodes the measured data signal y(l) by fitting a scanned signal 
model /(*, 0) to the measured data signal y{\) (shown generally as step 7) and then (in step 8) interpreting the encoded 
label data from the scanned signal model /(/, G). 

The scanned signal model y'(t, 0) is a model based on a typical or ideal measured label signal. The scanned signal 
model y(t, G) has a plurality of variable parameters 0 that varies the accuracy of the fit of the scanned signal model y(t, 
G) to the actual measured data signal y(t). Also, the accuracy with which the scanned signal model y(t, 6) fits the actual 
measured data signal y(t) depends on, e.g., the use of a fitting algorithm (shown generally as 7a). 

Once the scanned signal model y(f, 6) has been fitted appropriately to the actual measured data signal y{\), the 
final portion (step 8) of the decoding step 6 is to interpret the label data based on the values of the parameters ft used 
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in fitting the model to the scanned signal model j(t, e) to the measured data signal y(t). 

Referring now to FIGs. 2-3. a schematic diagram of a data reading system 10. e.g.. a barcode reading system, is 
shown. As discussed previously, an optical scanner (not shown) projects a bright spot 12 of light that moves across 
(shown generally as 14) a label 16 of interest, e.g.. a bar code label or other symbology data label such as printed 
characters, symbols or graphical information. A photodetector 22. with light collection lens 11, measures the time- 
varying intensity of backscattered light from label 16 and converts it to a corresponding time-varying measured or 
detected data signal y(t). e.g., an electrical signal current Then, as will be discussed in greater detail later, a bar 
code scanner or receiver 24 reads and processes the time-varying measured data signal y(t) to output a rectangular 
pulse train (shown generally as 26 in FIG. 2b) representative of label 16. From pulse train 26, the data elements or 
symbols encoded in the bar code label are determined in a conventional manner 

At a given point in time during reading of label 16 by photodetector 22, the measured data signal. >/(t). is the sum 
of elemental contributions of individual light intensities characterized by a point spread function h{l t r\, a, b, c) of light 
reflected from the point ft, x\) if beam spot 1 2 is centered at (a, b, 0) and the beam focal region is at "a distance c above 
label 16. Thus: 

>'(') = J| H^q^c) m(£,T|) dx\ 

A 

where a. b and care functions of time and where m(^ n ) is the printed source image. The point spread. function /?ft, 
y\, a, b, c) is assumed to be symmetrical and typically is modeled by a conventional Gaussian distribution density: 

25 (I'B) 2 •*• <n**>) 2 

hH.r\.a t b.c)=kB 2o(c)2 

where c(c) is proportional to the beam width and (a, b) are the coordinates of the beam center. By proper translation 
and scaling, and without loss of generality, the functions are assumed to be probability density functions. Thus, it is 
desirable to solve the inverse problem, i.e., to restore the undistorted image mft^i) from its distorted, scanned image 
e.g., measured data signal y(t). 

In the case of bar code data. assuming, e.g., constant velocity of beam movement across label 16 and a sufficiently 
large label 16 relative to beam spot 12, the measured data signal y(/) is modeled by receiver 24 as: 
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where s(x, X) is a general form of the input signal to be obtained from the measured data signal y{t), 6 = {a, X] is a 
parametric vector to be recovered from the measured data signal y(f), and 4>(f, c) is a Gaussian impulse response of 
backscattered light having the form: 



♦ (f.o) = 




where a is the point spread value, which is characteristic of the point spread function of the light spot 12 in the decoding 
so system 10. w y 

If a v .... a n is a binary representation of a bar code and X is a plurality of parameters including, e.g., a v a 2 , 
a n . .... t 1 and T, the general form input signal s(x, X) obtained from the measured data signal y(t) is defined as: 

55 1 

s(x 9 ^) = —2^a k p T {x-t k ), f4 = * + <* - 7)7 
1 k=i 
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where pj<x) is a rectangular pulse of width Tcorresponding to a bar in which pr(x) = 1 if Ixl £ 0.57, otherwise pj{x) : 
0. From these equations, it is seen that the measured data signal y(r) can be modeled linearly. 
Further derivation of the model j(t, 0) yields, generally: 



n 



4 = 1 



to where: 



75 



j 0.57 

w (r,7\a) = - {<|>(/ -jc,o*/x 



-0.57 



As shown generally in FIGs. 2-3, the scanned signal model /(/, 6) is used to generate pulse train 26 representative 
20 of label 16 from delected data signal y{() received by pholodelector 22 and processed by receiver 24. Such is accom- 
plished by determining coefficients a 7 , a^..., a n of detected data signal y[t) by fining scanned signal model /{(, 6) to 
detected data signal y{t). 

Conventionally, in fitting the scanner model y'(f, G) to the detected data signal y{t), signal processing methods of 
bar code reading systems assume that the relative blurring caused by light backscattered from adjacent bars of label 
zs 16 is insignificant. That is, a blur factor (a / 7) indicative of the relative point spread of light compared to the width of 
the label bars conventionally is assumed to be relatively small. In using such assumption, 0) a a^u(r -r^ 7, o) and 
conventional barcode recovery is accomplished, e.g., by sampling ;(r, 0) at moments t h or by using intersections of u 
(t - t k , 7, a) with its mean value, or by using zero crossings of its second derivative. 

However, intersymbol interference typically is not insignificant. Still, conventional bar code reading methods and 
30 systems fail to take its effects into account. That is, the contribution or effect of light backscattered from adjacent ■bar" 
or "space* elements a k _ ,u(t - t k . h T, o) and a k ^u (t - ( k ^ t 7, a) for /= 1, 2. ... on the overall decoding accuracy is not 
considered by conventional bar code reading systems. Additionally, other sources of interference and noise provide 
distortion to the measured data signal y(t). 

To fit the scanned signal model /(f, 0) to the measured data signal y(i) in the face of this blurring and other distortion, 
35 an optimization criterion is needed. An embodiment of the invention uses a deterministic expectation-maximization 
(DEM) algorithm to fit the scanned signal model /(*, 0) to the measured data signal y(t). 

The DEM algorithm is a novel generalization of the expectation-maximization (EM) algorithm, which is a conven- 
tional iterative method useful for finding maximum likelihood estimates of statistical model parameters from experi- 
mental data. Advantageously, the DEM algorithm is much more general and has more applicability than the EM algo- 
<*o rithm. 

Being more general than the conventional EM algorithm, the DEM algorithm provides solutions for what are known 
as positive inverse problems of finding parameters of non-linear models by fitting them toother models. Special cases 
of these problems include, e.g., approximation of a non-negative function with a parametric family of probability density 
functions, solving systems of non-linear equations, and solving systems of non-linear integral equations. 
45 Positive inverse problems are formulated, for example, in a manner of recovering an object s(x, X) from a parametric 

family of objects from its measured data signal y(t) having distorted image data, it the distortion can be modeled as: 



j{t$) = \s(X,\)h(x,t,\L)dx 



where h(x, t, u) is the point spread function, x\s an n dimensional vector, t is an m dimensional vector, and parameter 
0= (>-,u.) is also a multidimensional vector. In particular, if xand fare two-dimensional vectors (x, tGR 2 ), then the model 
j(l 0), which in data reading applications is referred to as the scanned signal model or the image distortion model, is 
55 interpreted as a distorted image of a planar object s(x, X). 

If parameter 0 is known, s(x, X) is capable ol being found as a solution of the previous integral equation. However, 
in the majority of practical cases, the parameter© is not known. Therefore, solutions are obtainable by minimizing the 
divergence between the measured data signal y(t) having distorted image data and its accompanying model ;(/; ft), e 
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g., by minimizing the Kullback-Leibler divergence: 



'(y|jb)=Jy(01og-^^f 

between the measured data signal y(t) and its model j(t, G). The Kullback-Leibler divergence is known to those skilled 
in the mathematical arts and is described, e.g., in the book S. Kullback, Information Theory and Statistics John Wilev 
& Sons, New York, 1959. y 
Without loss of generality, yfy and/ft G) are assumed to be probability density functions. Since yf/; does not depend 
on 6, minimization of K(yUjQ) is equivalent to the maximization of the mean log-likelihood function: 

L MU)= Jy(01og j{t3)dt . 

The solution of the inverse problem is determined by 6 = arpmax L(y ll/ 8 ). Although this equation typically is difficult 
to solve directly, the DEM algorithm is useful for finding the maximum of the function iteratively. More specifically, the 
DEM algorithm is based on presenting L(B) in the decomposed form L(G) = O(0, e p ) - H(e. 0 p ), where O(0. 0 ) and H 
(0, G p ) are auxiliary functions of the observed data y. If G^ f can be found such that 0(G p+r e p ) > O(0 , 0 ) and H(Q 
G p ) < H(0 p . G p ),then L(6 p+T )> L(G p ). In this manner, a sequence 9,, 6^,... G p is created that increases L(Q) monotonicaNy 
The sequence converges to a maximum of L(G) if the functions and their defined region satisfy very mild regularity 
conditions, which are known to those skilled in the art and are found, e.g., in the book W. I. Zangwill, Nonlinear Pro- 
gramming: A Unified Approach, Prentice Hall, Englewood Cliffs, NJ, 1969, pp. 91-92. 

H "(6 p , G p ) is a global maximum of H(G, G p ), then L(Q) increases only if O<0, 0 p ) increases. Thus, the DEM algorithm: 

monotonically increases Z_(G) and finds its maximum under mild regularity conditions 

The auxiliary functions O(0, G p ) and H(G, G p ) are not uniquely identified. Their selection is based on the simplicity 
of finding a maximum of 0(G, G p ), which directly affects the algorithm convergence rate. The auxiliary functions is found 

e.g., by: 

!°9 A*> G) = log q(z> G) - log p(z, G) 

where z = (s, x) is an augmented variable, p(z, G) is its probability density, and q(z, G) = p(z, G) /(/, 0). Then the auxiliary 
functions are constructed as: ' 



<2(e ,e p ) = J J iog{q( z ,Q ) } p(z,e p ) y (x)ds 



Since H(9. e p ) attains its global maximum at e = e p , the DEM algorithm is applied to find a maximum of L(jMJ e ) iteratively 
Alternately, since a maximum of L(0) is a minimum of -L(Q), the DEM algorithm is uselul tor finding a minimum 

Therefore, the DEM algorithm is adaptable to applications where a minimum rather than a maximum is sought 

The conventional EM algorithm is one specific application of the DEM algorithm. Specifically, when there are ob- 
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servations x v x 2 . ... x N oi a random, mdimensional variable x with the probability density function /(/, 6), the EM algorithm 
is obtained from the DEM algorithm when y(l) is defined as: 

N 



k=\ 



70 where 6(/ - x^) is the Dirac delta-function, which is known in mathematical physics and applied mathematics. 

When the DEM algorithm is used in, e.g., bar code reading methods and systems, the auxiliary function 0(8, e p ) 
is expressed, e.g., as: 



G(e,e p )= J jiog{^(t-x J a)s(x))w(t,x 9 e p )dtdx 

where w{t, x, e p ) = <|>(f - x, o) s p (x) y[f) / 6 p ). 

II the point spread value o is known, the scanned signal model /(t, 0) is a linear function of s(x) and it is not difficult 
for those skilled in the art to find a maximum of function 0(6, 6 p ). The result is the following EM algorithm (i.e., a specific, 
"linear" version of the inventive OEM algorithm): 

■V.U)= jw(/^ i e / ,)rfr = 5 / ,(x)jr(r,6 p )<t)(/-^o)^ 

30 where r[t, fl) = y(0 //(', G p ). This particular algorithm was introduced by Vardi and Lee (e.g., see Y Vardi and D. Lee, 
"From image deblurring to optimal investments: maximum likelihood solutions for positive linear inverse problems," 
Journ. Royal Statist: Soc, Vol. 55; 1993, pp. 569-612) for solving general linear inverse problems (i.e., problems un- 
related to recovering bar code symbols and other symbology data). 

FIG. 4 illustrates the results of an application of the Vardi and Lee EM algorithm solving a linear inverse problem 

35 (using a known point spread value o) after approximately TOO iterations. As shown, the transitions of the measured or 
received data signal y(t) 32 (solid line) is moderately well modeled by a scanner model signal 34 (dashed line) fitted 
thereto. Furthermore, scanner model signal 34 depicts what would be a deblurred measured signal and outputs a 
partial rectangular pulse train useful by subsequent processing in decoding particular symbology characters therefrom. 
However, the Vardi and Lee "linear" version of the algorithm provides no useful solution to non-linear inverse 

*o problems. That is, the Vardi and Lee algorithm is not helpful when the point spread value a is unknown and, as mentioned 
previously, the point spread value a usually is not known, because it depends on the scanner distance from the bars 
in the bar code label, which distance is a random variable. 

When the point spread value a is unknown, the problem becomes non-linear and, conventionally, a unique solution 
for s(x) is not available. However, with the use of a data-characteristic model, e.g., ;(/, 6), the invention obtains a unique 

45 solution for non-linear inverse problems by fitting the data-characteristic model to the collected data using the DEM 
algorithm and determining the data, i.e., solving fors(x), based on the parameter values of the data-characteristic model. 

Accordingly, for applications such as decoding bar code label data and other symbology data, the use of the 
scanned signal model /(t, 6) enables unique solutions to be obtained for the collected data. Although a still is unknown 
and the inverse problem remains non-linear, the number of parameters needed to express the data-characteristic 

so model is limited sufficiently so that unique solutions are available for the collected data when the data-characteristic 
model is fitted thereto. 

Thus, a form of the DEM algorithm more general than the EM algorithm is useful for solving non-finear inverse 
problems. Specifically, the DEM algorithm is obtained by maximizing the function O(0, 6 p ), where the parameters a, o, 



t and Tare defined as: 
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a 2 ^,= J j(t-x) 2 w(t,xj p )dtdx 



{'..,*. • ^, } = arg max £ log] \ J w(r, *,6 , ) 

* = l I * J/, -0.57 

20 Typically, / 7<p ^ 7 and T p ^ 1 are difficult to determine from these equations. However, by using an approximation u{( -/^ 
T, a) as 0(r -/^ a) and the relationship t k = f 7 + (/c - 7)T the inventive DEM algorithm is obtained by maximizing the 
function 0(6, 6 p ), where the parameters a 7 , a 2 ,... ( a n , o, f 7 and Tare defined as: 



25 a *.P+1 ~ a k.p V k,p 



where 

40 
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In an alternative embodiment of the invention, the use of the DEM algorithm as part of the bar code recovery 
method and system as described herein is supplemented by the use of conventional zero crossings determinations, 
which typically are used alone in some conventional bar code recovery methods. Specifically, the conventional zero 

10 crossings technique is used to determine the initial values (i.e., for p = 0) of the parameters G 0 = \t 10 , T 0 , o 0 , a 1 0 , 
a 2 Q, ■ a n0 ). N these parameter values are suitable for purposes of fitting the scanned signal model /(/, 0) to the 
received data signal y(t) (i.e., if the recovered bar code from the scanned signal model fit is determined to be a proper 
bar code), then no further signal processing is needed. Otherwise, the initial zero crossings parameters 0 O = { t 1 0 , T& 
<3 0 , a 1 0 , a 2 o a n.o) are used on, y as ' n ' l 'al parameter estimates for the DEM algorithm. 

J 5 For example, FIG. 5 depicts a received data signal y(t) 42 as compared to a corresponding scanner signal model 

fitted signal 44 that has used parameter values determined by a conventional zero crossings technique. In this specific 
case : a recovered bar code based on the scanner signal model fitted signal 44, as shown, would be deemed improper. 

By comparison, FIG. 6 depicts a scanner signal model fitted signal 54 (after 40 iterations) to a received data signal 
y(t) 52 using the DEM algorithm. Moreover, as shown in FIG. 7, when the zero crossings parameters e 0 = {t 1 q, o 0 , 

20 a j q, a 2 o . 3 n qI are used as initial parameter estimates for the DEM algorithm, a scanner signal mode! fitted signal 
64 satisfactorily models a received data signal y(t) 62 after only 11 iterations. In such case, the use of the DEM algorithm 
in fitting the scanner signal model j(t, 9) to the received data signal y(t) is improved (in terms of both accuracy and 
number of iterations required) when parameters determined by a zero crossings technique are used as the initial 
parameter estimates for the DEM algorithm. 

25 it is to be understood that the DEM algorithm has many applications not discussed specifically herein but that are 

obvious to those skilled in the art. For example, the DEM algorithm is useful in the minimization of the number of states 
of hidden Markov models used in speech and handwriting recognition., the optimization of the structure of natural 
language models, and for communication channel bljnd equalization. 

Thus, post -computational activity can be of a variety of types, but will generally be similar to corresponding prior 

30 art activity. Exemplarily, it includes, for example, the display of information (e.g., a price of the item having a bar code 
label) or a change in the state of a physical system (e g., opening a security gate). 

It will be apparent to those skilled in the art that other changes and substitutions can be made to the symbology 
data decoding method and system described herein without departing from the spirit and scope of the invention as 
defined by the appended claims and their full scope of equivalents. 

35 

Claims 

1. A method for determining characteristic information from a signal, said method comprising the steps of: 

40 

providing (5) a measured data signal y(t) representative of a plurality of input data capable of including distor- 
tions therein; 

fitting (7) a scanned signal model /(/, 0) having a plurality of parameters 0 to said measured data signal y(t) 
by adjusting the values of at least one of said parameters 0, said scanned signal model y(t, 0) capable of having 
4 $ a nonlinear relationship to said plurality of parameters 0, 

said fining step comprising using a DEM algorithm to determine the values of at least one of said plurality of 
parameters 0; and 

determining (8) said symbology information from said scanned signal model /(/, 0) based on the values of said 
parameters 0 when said scanned signal model j(X, 0) is fitted to said measured data signal y(t). 



so 



55 



2. The method as recited in claim 1 , wherein said DEM algorithm is: 

0 p+J = arg max O(0, 0 p ), for p = 0, 1 

where O(0, 0 p ) is an auxiliary function of said parameters 0. 

3. The method as recited in claim 2. wherein said fitting step uses said DEM algorithm to obtain an input signal s(x) 
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from said measured data signal y(t), wherein: 



(2(6,6,)= J jlog{<W- x t O)s{x))w(t,x,e P )dtdx , 

where w(t, x, fl p ) = W - x, o) s p (x) y(t) /j(t, e p ), and where W. o) is a Gaussian impulse response of backscatlered 
light having the form: 

10 



1 2o 2 



75 

where o is a point spread value. 



4. The method as recited in claim l , wherein said fining step uses said DEM algorithm to obtain an input signal six 
a) defined as: ? v ' 



20 



25 



40 



45 



50 
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s(x,\) = —^a k p T (x-t k ), 

from said measured data signal y(Q. said input signal s(x. X) having a plurality of parameters capable of being 
solved by said DEM algorithm, said plurality of parameters defined as: . 

30 a *.P*i = a k.P V k. P i 

n 

35 * = I 



T P*i = ( M 1. P - ApWiBp- A?) 



where 



A p = X (* - 1K, +1 , B p = 2 (A: - 1)« 

*=2 A=2 



* = 2 



10 
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v ..p = ] < x - '«>- c > 1 p^ dx 



C * p+l = — ?(*- t pi f y(xM (* - t ip , a ,, ) / j(x-e p )dx. 

10 V- 

5. The method as recited in claim 1, wherein said fitting step includes minimizing the Kullback-Leibler divergence 
defined as: 

75 

y(0 



20 



K(y\\j 6 )= jKOlog-^// 



6. The method as recited in claim 1, wherein said titling step includes maximizing a mean log-likelihood function 
defined as: 

25 



Uy\\j B )= ]y(t)logj(t,Q)dt 



30 



7. The method as recited in claim 1 , wherein said fitting step further comprises using zero crossings of said measured 
data signal y[\) to determine initial values for at least one of said plurality of parameters 6, said initial parameter 
values usable by said DEM algorithm as initial estimates for said parameter values. 

35 

8. The method as recited in claim 1 , wherein said determining step further comprises generating a pulse train corre- 
sponding to said scanned signal model 6) when fitted to said measured data signal y(t) and identifying said 
characteristic information based on said pulse train. 

40 9. The method as recited in claim 1 , wherein said scanned signal model is defined as: 



45 



j(t,Q) = \<b(t-x,o)s{x,X)dx 



where s{x, X) is an input signal obtainable from said measured data signal y(t), e = {o.X] is a parametric vector 
recoverable from said measured data signal y(t), and where 4>(f, o) is a Gaussian impulse response of backscattered 
50 light having the form: 

^(/,c) = - 7 ^e 2 ° 2 
55 j2no 

where ct is a point spread value. 
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10. The method as recited in claim 1. wherein said scanner model is defined as: 



jV,Q) = £a t u(r-/*,7\a), 

where: 



j O.ST 

u(r r 7\a) = - j<p(t-x,o)d. 



-GST 

and where <!>(/, a) is a Gaussian impulse response of backscattered light having the form: 



J2ncs 



and where o is a point spread value. 



11. A system for reading bar code labels, said bar code labels having at least one data character encoded therein, 
said bar code label said system comprising: 

a photodetector (22) for collecting light backscattered from a bar code= label, said backscattered light corre- 
sponding to said bar code label; and 

a receiver (24) in operable connection with said photodetector for converting said backscattered light into a 
time-varying measured data signal y(t), fitting a scanned signal model y(/, 9) having a plurality of parameters 
G to said measured data signal y(t) using the DEM algorithm, and determining said data characters encoded 
in said bar code label from said scanned signal model ;(/, 9) based on the values of said parameters 9 when 
said scanned signal model /'(/, 9) is fitted to said measured data signal y(1). 
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(54) Method and system for decoding distorted image and symbology data 



(57) Disclosed is a system and method for reading 
and decoding distorted data, e.g. . bar code label or other 
symbology data, by fitting a scanned signal model j(U 
8), using a Deterministic Expectation-Maximization 
(DEM) algorithm, to the measured or detected data sig- 
nal y[t). The DEM algorithm, which is a novel variation 
of the conventional Expectation-Maximization (EM) al- 
gorithm, enables the system and method according to 
the invention to determine the values of a plurality of 
parameters 0 defining the scanned signal model /(/, 6) 
as it is being fitted to the measured data signal y(Q. 



Based on these values, characteristic information such 
as data elements or symbols encoded within the image 
data are generated using conventional signal process- 
ing techniques. The use of the DEM algorithm advanta- 
geously deblurs the signal by reducing the effects of in- 
tersymbol interference. In an alternative embodiment, 
the DEM algorithm is supplemented by a conventional 
zero crossings decoding technique. In this manner, the 
zero crossings technique is used to obtain initial param- 
eter estimates for the DEM algorithm in determining the 
parameters 8 of the scanned signal model j[t, 6). 
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